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ABSTRACT 

We have re-analyzed archival HST and Subaru data of the recently discov- 
ered planetary mass companion (PMC) GQ Lup B. With these we produce 
the first R and I band photometry of the companion and fit a radius and ef- 
fective temperature using detailed model atmospheres. We find an effective 
temperature of 2335 ± 100K, and a radius of 0.38 ± O.O5R and luminosity 
of log(L/L ) = —2.42 ± 0.07 (at 140pc). Since we fit wavelengths that span 
most of the emitted radiation from GQ Lup this luminosity estimate is robust, 
with uncertainty dominated by the distance uncertainty (±50pc). The radius 
obtained for 140pc (0.38R©) is significantly larger than the one originally de- 
rived and larger than model predictions. The mass of the object is much more 
model-dependent than the radiative properties, but for the GAIA dusty models 
we find a mass between 10-20 Mj up , in the range of the brown dwarf and PMC 
deuterium burning boundary. Assuming a distance of 140pc, observations fit 
to la the Baraffe evolution model for a ~ 15 Mj up brown dwarf. Additionally, 
the F606W photometric band is significantly overluminous compared to model 
predictions and other brown dwarfs. Such overluminosity could be explained by 
a bright Ha emission from chromospheric activity, interaction with another un- 
detected companion, or accretion. Assuming that GQ Lup B has a bright Ho; 
emission line, its Ha emission strength is 10 _1 ' 71±0 ' 10 Lb o i, significantly larger than 
field late-type dwarfs. GQ Lup B might be strongly accreting and still be in its 
formation phase. 

Subject headings: stars: imaging, stars: pre-main sequence, stars: low-mass, 
brown dwarfs, (stars:) planetary systems, techniques: photometric 
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Introduction 



Direct exoplanet detection around stars is a challenging endeavor, but po ssible for com- 



panions to the youngest stars ( Neuhauser et al. 2005 ; Chauvin et al. 2005al lbh . Unlike the 



exoplanets found by precision radial velocity techniques that provide little information (albeit 
model independent) about the physical properties of planetary mass companions (PMCs), 
the physical properties of these PMCs must be inferred by comparing atmosphere and evo- 
lution models to observed spectra or photometry. 



One of the most recent PMC candidates identified orbits the star GQ Lu p (INeuhauser et al 



2005) (N05). This star is a young K7eV TTauri star in the Lupus I cloud (jTachihara et al. 



19961 ). The star has an estimated age of less than 2Myr (N05 ) and is situated at 140pc with 
a potential range fr om 90 to 190pc ( Wichmann et al. 19981 ; Neuhauser fc Brandner 1998 



Knude fc Hoall998l ) . The original analysis of the GQ Lup PMC candidate, called GQ Lup 
B, is based on K- and L'-band photometry along with a K-band spectrum. From this early 
work, an effective temperature of 2050 ± 450K, a radius of 0.12 ± 0.06R m , a luminosity of 
loz(L/Lr*)__= - 2.37 ± 0.41 and hgg of 2.52 ± 0.77 are obtained JNeuhauser et all 12005 



Neuhauser! 120051 1 . 



Since other wavelength bands are available from HST and Subaru that span most of 
the radiation emitted by GQ Lup B, it is possible to conduct a more complete photometric 
analysis. This analysis is presented below along with the implications of a detected R-band 
overluminosity. 



2. GQ Lup Photometry 



In addition to the (N05) VLT data, GQ Lup B has been previously observed by both 
the Subaru telescope, program o02312, and the Hubble Space Telescope (HST), programs 
SNAP7387 and 9845. The coronagraphic imager with adaptive optics instrument (CIAO) 
(IMurakawa et al.l 120041 ) was used at the Subaru telescope while the Wide Field Planetary 
Camera No. 2 (WFPC2) and NICMOS were used with HST. 



2.1. HST Visible and NIR Photometry 

Data were retrieved from the public MAST Hubble Space Telescope archive at STScI 
for the filters F606W, F814W, F171M, F190N and F215N using the automated reduction 
pipeline. The companion is clearly visible 0.7" West of GQ Lup in all filters. The GQ Lup 
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PSF in each filter was first subtracted using reference PSFs of a second star ob served in th e 
same program or using simulated PSFs produced by the Tiny Tim software (IKristl Il993l ). 
selecting the one that gives the smallest residual at 0.7" separation. The simulated PSFs in 
each filter were then used to estimate the GQ Lup B flux. PSFs were simulated with five 
times the sampling, shifted, and binned to the detector resolution. The set of parameters, 
i.e. fractional pixel PSF position and flux normalization, that minimize the RMS noise inside 
a 6 x 6 pixels box centered on GQ Lup B was kept. Magnitude errors were estimated by 
calculating the RMS value of the total GQ Lup subtracted PSF residual flux at the same 
angular separation using the same 6x6 pixels box but at a different field angle. Regions 
contaminated by residual flux from the diffraction spider were avoided. Total companion 
fluxes were obtained by integrating the Tiny Tim simulated PSFs that b est subtract the 



GQ L up B PSF. The WFPC2 charge- Transfer efficiency bias was corrected (jWhitmore et al. 



19991 ). but the amplitude of the effect is small, less than 5%, since GQ Lup is bright (more 



than 3000 counts inside a 2 pixel radius aperture). Table [T] shows the obtained apparent 
magnitudes and est imated error bars. Inters tellar extinction, assumi ng A v = 0.4 ± 0.2 
(IBatalha et al.l 120011 ) and the extinction law of iRieke fc Lebofskyi (119851 ). is also tabulated. 



2.2. Subaru Photometry 

Using the public Subaru archive (SMOKA), data for K, CH 4 and L' filters were re- 
trieved. For K and CH 4 , since a dithered pattern was used to acquire GQ Lup, a sky image 
was constructed from the median of all acquired images in each band. For L', the sky image 
sequence was used to subtract the thermal background. Since no flat field images are avail- 
able, images were simply registered at the image center, median combined, and a 180 degrees 
rotation was used to subtract the smooth PSF halo and residual sky background (an angle 
of 170 degrees was used for CH 4 data due to a bad pixel located at 180 degrees of GQ Lup B 
position). GQ Lup B is clearly visible in all three filters and was not saturated or occulted 
behind a coronagraph. The PSF average peak flux per readout in raw images for the CH 4 , K 
and L' bandpasses are respectively 30,000, 40,000 and 50,000 electrons. PSF peaks are thus 
below the 60,000 electrons level (CIAO Aladdin-2 detector has a 5% non-linearity at 60,000 
electrons) in all bandpass images, detector non-linearity is thus neglected. Therefore, Amag 
measurements were possible in each filter. GQ Lup B magnitude differences were found by 
optimizing the PSF subtraction from the PSF of GQ Lup. Magnitude errors were estimated 
by calculating the flux RMS variation of 10 square boxes of 1.5 X/D width at 10 different 
angular positions and at the same separation as GQ Lup B. Apparent magnitudes were de- 
termined using known H (to determine the CH 4 magnitude), K and L' m agnitudes for GQ 



Lup, H=7.70±0.03 mag, K=7.10±0.02 mag (2MASS) and L=6.05±0.13 flGlass & Penston 
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19741 ; iHughes et al.lll994l ). Magnitude differences between the L' and the L filters are negli- 



gible. Our L' magnitude differs by 0.74 magnitude compared to the one published by N05 
while the K (Ks for N05) magnitudes are consistent to 1 a accuracy (see Tabled]). A slice 
of the L' image through GQ Lup A and B is given in Fig. [H the N05 measurement is likely 
contaminated by flux from the primary. 



3. GQ Lup B Temperature and Radius Fit 

A solar metalicity subset of the GAIA dusty model atmosphere grid (Hauschildt et al., 
in preparation) was used to fit the photometry in Sect. [2] and determine the radius and 
temperature of GQ Lup B. The fact the all the data was not acquired simultaneously but 
spread over many years could suggest a problem related to a potential GQ Lub B (and A 
for contrast ratio measurements) variability. Since the HST, VLT and Subaru H and K 
photometry were acquired two years apart and agree to ~ la suggest that GQ Lup B (and 
A) variability is small for our dataset. Given that surface gravity has only a small effect 
on predicted broad band photometry and that such young low-mass objects typically have 
low gravities, a standard logg = 3 was assumed. This gravity is consistent with the one 
found by N05 from a low resolution K-band spectrum, i.e logg = 2.52 ± 0.77. The model 
fluxes were convolved and integrated over the appropriate transmission curves for each filter 
listed in Table 1. For the F814W HST filter, the total optical system transmission was used 
instead of the filter transmission due a drop of sensitivity in a spectral region where GQ Lup 
B shows a significant increase of luminosity. These synthetic broad-band flux densities were 
interpolated to produce a uniform square grid in radius and temperature space. Observed 
magnitudes were corrected for distance and estimated interstellar reddening and transformed 
to fluxes using a calibrated Vega spectrum. The shape of the observed broad-band SED from 
R to L-band determines the effective temperature while the overall scaling needed to best 
match the observed fluxes (for a fix distance) determines the radius. A standard deviation 
for each model, cx mo dei, was calculated using the following standard equation, 



where N is the number of photometric measurements and o, e and m are respectively the 
observe data point, its error bar, and expected model-derived flux. The F606W bandpass was 
not included since this point is clearly overluminous compared to the near-infrared bands (see 
Sect. HJ). Figure [2] shows the temperature as a function of radius for GQ Lup B. Error contours 



^modcl 



\ 



- 5 - 



are defined as models that de viate by 1, 3 and 5 a from the best fit model. The predicted 
temperature and radius of the iBaraffe et al.l (120031 ) evolution models for 1 and 5 Myrs and 7, 
10, 12, 15 and 20 Mj up are also shown. The best fit found is at 0.65cx from the model with an 
effective temperature of 2335K and radius of O.38R (assuming a distance of 140pc). The fit 
was then rerun for both maximum and minimum allowed distance and interstellar extinction 
to estimate error bars (see Table [2]). The overall fit to account for both the distance and 
interstellar extinction is thus an effective temperature of 2335±100K and a radius of range of 
0.38±0.05<ii4o R.0, where <ii 40 = d/140 and d is the distance in parsecs (see Fig. [3]). Note that 
the fit would be equally consistent with a log g = 4 in agreement with the fact that broadband 
photometry is insensitive to logg (the mass of GQ Lup B can thus not be contrained by our 
adopted logg and derived radius). The resulting luminosity, accounting for the correlated 
error bars, is equal to log(L/L Q ) = —2.42 ± 0.07 + log(df 40 ). Interstellar extinction does 
not significantly change the derived effective temperature but the distance adds a possible 
range three times the derived radius error bar. A better distance estimate would be required 
to further constrain the radius and the luminosity of GQ Lup B. If the GQ Lup system is 
really at 140pc, for our best temperature fit, our radius measurement is larger (by 3a) than 
model predictions (0 .32 R ), a result similar with what is observed in young star clusters 
(iMohanty et al.ll2004l ). The larger than expected radius could also be explained if GQ Lup B 
is actually an unresolved binary. At a la confidence level, our derived effective temperature is 
consistent but has an error bar five times smaller than the one derived by N05 and our radius 
is more than three times larger. Note that the radius (0.12 ± 0.06R P ), effective temperatur e 
(2050 ±450) and luminosity (Ior( L/L m ) = -2.37 ± 0.41) given in iNeuhauser et all koO$ \: 
Neuhauserl (120051 ) ; lGuenther et al.l (120051 ) are incompatible; for such radius and temperature, 
the luminosity should be closer to log(L/L ) = —3.67, clearly inconsistent with the observed 
photometry (see Fig. [3]). It is clear from Fig. [3] that the F606W magnitude is overluminous 
compared to model predictions - the implications of this are discussed below. GQ Lup B is 
also overluminous in the R-band when comp ared to other brow n dwarfs (see, for example, 
the visible spectra of Pleiades brown dwarfs ( Martin et al. 200ol )). 



As discussed in N05, planetary evolution models are hig hly uncertain f o r you ng objects 



However, it is still interesting to compare our results to the IBaraffe et al.l (120031 ) evolution 

~lMyr objects 



(Baraffe et al. 


2002) 


Mohantv et al 


2004 



(Baraff e et al.ll2002r). Such problems have been confirmed by o bservations (IKonig et al.ll2002 



range (90 to 140pc) and age (< 2 Myr), our temperature and radius are consistent with a 
10-20 Mj up PMC/brown dwarf at la, at the boundary bet ween a PMC and brown dwarf 



similar to the AB Pic PMC candidate (IChauvin et al.ll2005bl ). Assuming a distance of 140pc 



a slightly higher effective temperature of 2400K and a smaller radius, our observations are 
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consistent with a ~ 15 Mj up brown dwarf (Baraffe models) to la accuracy. 



Discussion 



The HST F606W magnitude is significantly overluminous compared to model predictions 
(three magnitudes, see Fig [3]). Such overluminosity can be explained by an unmodeled effect 
such as a bright Ha emission line from chromospheric activity, interaction with another 
undetected companion, flaring, accretion or some other unknown process. Assuming that the 
observed F606W overluminosity is coming from Ha emission, we can estimate the strength of 
emission by calculating the log ratio of its Ha luminosity to its bolometric luminosity. Using 
the simulated flux normalized spectrum for GQ Lup B, the F606W filter bandpass profile 
and assuming that all the observed flux in the F606W filter comes from a bright Ha emission 
line, we find GQ Lup B Ha emission strength log(L#a,/Lboi) of —1.71 ±0.10. If this emission 
if from ac cretion, it is s t rong enough to be classified as accreting following the criterion 
defined in iBarrado et al.l (120031 ) . This observed emission strength is signi ficantly larger, by 
an order of magnitude, than what is found for field M, L and T dwarfs (|Gizis et al. 2000 ) 



or even in peculiar late type dwarfs (ILiebert et al.l 119991 ; iBurgasser et d.ll2000tl2002n. Such 
pecul iar dwarfs are thought to be young low mass objects, ~10 Myr 3-20 M Jup (ILiebert et al. 



20031 ). though not as young as GQ Lup ; GQ Lup B could be a very young example, still 
bound with its primary, of such objects. 

A visible spectrum of GQ Lup B is needed to confirm the Ha emission. If the Ha 
emission is present, the emission li ne 10% width can be used to discrimina te between ac- 
cretion and chromospheric activity (IMuzerolle et al.ll2003l ; iNatta et al.l 12004 ). Detection of 
Pa/3 and B ry lines in the near-infra r ed could also be used t o confirm accretion. Br7 was not 
detected in lNeuhauser et al.l (120051 ); iGuenther et al.l (120051 ) K-band spectrum, but since this 
line is harder to dete ct and fainter tha n the Pa/? line, such non-detection does not imply 
no ongoing accretion (INatta et al.l 120041 ) . A time series photometric/spectroscopic analysis 
could distinguish between transient emission due to a strong flare or constant emission char- 
acteristic of accretion. Only two R-band images acquired less than four minutes apart are 
available from the HST archive. The individual photometry of these two images agrees to 
a few percent, consistent with the constant emission from the accretion hypothesis. Due to 
the limited time sampling of these observations, more data with bigger time intervals are 
needed to better quantify the Ha flux variability. Searching for an eclipse could also confirm 
the interacting binary hyp othesis. Strong Ha e mission could also be a sign of a runaway 
accretion as postulated by iFortney et al.l (120051 ) for the core accretion-gas capture model, 
although such observations would be very fortuitous due to the relatively short timescale ex- 
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pected for this phase. Another more probable possibility is that GQ Lup B simply has form 
through independent collapse of a molecular cloud fragment at its current separation and is 
still accreting its circumstellar disk. A PMC search at smaller separations could distinguish 
the two possibilities. If GQ Lup B is still in its formation phase, it would become a unique 
object to study ongoing PMC formation. 

If the Ha emission is from accretion, an unresolved/undetected circumstellar disk must 
be present around GQ Lup B. Our photometric analysis did not detect an infrared excess at 
L', although the L' photometric point is ~ la overluminous compared to the models. 



5. Conclusion 



We have reanalyzed available Subaru and HST data to fit GQ Lup B radius and effective 
temperature using model spectra. Our derived effective temperature (2335±100K) is slightly 
hotter than the one derived by N05 but with a substantially smaller error bar. Our derived 
radius for 1 40pc (0.38 ± 0.05R Q ) is more than three times larger than what was found by 
Neuhauserl (120051 ) and larger than model predictions. At that distance, our result is consistent 
to la accuracy with a ~15 Mj up brown dwarfs (Baraffe models). If GQ Lup B is confirmed 
to be strongly accreting, it might be a young, still forming/ contracting PMC or brown dwarf. 
A better distance estimate for the GQ Lup system is needed to further constrain GQ Lup B 
radius and its luminosity. At the one sigma confidence level, our measurement is consistent 
with a 10-20 Mj up companion, at the mass boundary between PMC and brow n dwarf. Our 
results are in agreement with the independent work of iMcElwain et al.l (120071 ). 



Based in part on data collected at Subaru Telescope, NAO of Japan, and with the 
NASA/ESA HST, obtained from the STScI data archive (operated by the AURA, NAS 5- 
26555). This publication makes use of data products from 2MASS, which is a joint project of 
the University of Massachusetts and the IPAC/Caltech, funded by the NASA and NSF. This 
research was performed under the auspices of the DOE by the UC, LLNL under contract 
W-7405-ENG-48, and also supported in part by the NSF Science and Technology CFAO, 
managed by UCSC under cooperative agreement AST 98-76783. The authors thank Eric 
Becklin & Ben Zuckerman for helpful discussions. 
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Table 1: GQ Lup B Photometry 



Filter 


Date acquired 


Vega 


mag 


mag (Neu05) 


Interstellar 






zero point 






Extinction 


F606W 


1999-04-10 


22.92 


19.19 ±0.07 




0.4 ±0.2 


F814W 


1999-04-10 


21.67 


17.67 ±0.05 




0.2 ±0.1 


CH 4 off 


2002-07-17 




13.76 ±0.26 




0.07 ±0.04 


F171M 


2004-04-20 


20.19 


13.84 ±0.13 




0.07 ±0.04 


F190N 


2004-04-20 


18.48 


14.08 ±0.20 




0.06 ±0.03 


F215N 


2004-04-20 


18.25 


13.40 ±0.15 




0.05 ±0.02 


K 


2002-07-17 




13.37 ±0.12 


13.1 ±0.1 


0.05 ±0.02 


V 


2002-07-17 




12.44 ±0.22 


11.7 ±0.3 


0.02 ±0.01 



- 11 - 



Table 2: Best Radius and Temperature Fit 





Dist. (pc) 


IEy 


Teff(K) 


R(R ) 




90 


0.4 


2335 ± 80 


0.25 ±0.03 


Dist. 


140 


0.4 


2335 ± 80 


0.38 ±0.05 




190 


0.4 


2335 ± 80 


0.52 ±0.07 




140 


0.2 


2320 ± 80 


0.38 ±0.05 


IE 


140 


0.4 


2335 ± 80 


0.38 ±0.05 




140 


0.6 


2360 ± 80 


0.38 ±0.05 


Final 


140 


0.4 


2335 ± 100 


0.38 ± 0.05cii4o 
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Fig. 1. — GQ Lup B L' magnitude measurement. Solid line shows GQ Lup PSF intensity 
profile going through GQ Lup B. Dotted vertical line shows the separation of GQ Lup B. 
Dashed line shows the same intensity profile after sky subtraction and GQ Lup PSF halo 
subtraction. The two horizonthal dotted lines show N05 and our magnitude estimates. 
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Fig. 2. — Temperature-radius fit for 140 (left) and 90 and 190pc (right) distance. The 
diamond symbol represents the best fit. Three contour levels at 1, 3 and 5a from the best 
fit are shown. Model predictions for 7, 10, 12, 15 and 20 Mj up and 1 (dotted line) and 5 Myr 
(dashed line) are also shown. 
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Fig. 3. — Best temperature-radius fit (solid line) for 140pc, logg = 3, typical 0.4 IE in 
V band and 1 Myr. Luminosity is log (L/L ) = —2.42. Dashed line is the same fit with 
log# = 4. 



